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ABSTRACT 

We have detected the S(l), S(2), S(3), S(4) and S(5) pure rotational lines of 
molecular hydrogen toward the outflow source HH 54, using the Short Wavelength 
Spectrometer (SWS) on board the Infrared Space Observatory (ISO). The observed 
H 2 line ratios indicate the presence of warm molecular gas with an H 2 density of at 
least 10 5 cm -3 and a temperature ~ 650 K in which the ortho- to para-H 2 ratio is 
only 1.2 ± 0.4, significantly smaller than the equilibrium ratio of 3 expected in gas at 
that temperature. These observations imply that the measured ortho- to para-H 2 ratio 
is the legacy of an earlier stage in the thermal history of the gas when the gas had 
reached equilibrium at a temperature ,$ 90 K. Based upon the expected timescale for 
equilibration, we argue that the non-equilibrium ortho- to para-H 2 ratio observed in 
HH 54 serves as a chronometer that places a conservative upper limit of ~ 5000 yr on 
the period for which the emitting gas has been warm. The S(2)/S(l) and S(3)/S(l) 
H 2 line ratios measured toward HH 54 are consistent with recent theoretical models of 
Timmermann for the conversion of para- to ortho-H 2 behind slow, ‘C’-type shocks, but 
only if the preshock ortho- to para-H 2 ratio was 0.2. 


- 2 - 


1. Introduction 

Molecular hydrogen is the dominant constituent of dense interstellar gas, and Ho quadrupole 
transitions are a significant coolant of such gas over a wide range of temperatures and densities 
(Neufeld, Lepp & Melnick 1995). Although quadrupole vibrational transitions of Ho near 2 /tm 
have been widely observed from interstellar clouds that are strongly irradiated by ultraviolet 
radiation or have been heated by shock waves to temperatures T ~ 2000 K (e.g. Shull & Beckwith 
1982; Beckwith et al. 1978; Gatley et al. 1987; Chrysostoniou et al. 1993; and many others), 
molecular hydrogen within cooler clouds that are not UV-irradiated is only detectable by means 
of its pure rotational lines. With a broad spectral coverage that is uninterrupted by atmospheric 
absorption regions, the Short Wavelength Spectrometer (SWS; de Graauw et al. 1996) on board 
the Infrared Space Observatory (ISO; Kessler et al. 1996) has allowed the complete H 2 rotational 
spectrum to be observed from the interstellar medium for the first time. 

Pure rotational emissions from H 2 have been detected by ISO toward several regions of star 
formation within the Galaxy - including S140 (Timmermann et al. 1996), Cepheus A West (Wright 
et al. 1996), the BD+40°4124 group (Wesselius et al. 1996), and Orion-IRc2 (van Dishoeck et 
al. 1998) - as well as towards several external galaxies (e.g. Arp 220; Sturm et al. 1996). The 
observed H 2 line ratios provide valuable constraints upon the physical conditions in the emitting 
gas, implying typical temperatures of 500 - 800 K for those Galactic sources that have been 
observed to date. These temperatures are much greater than those typical of cold quiescent clouds 
(~ 30 K), and presumably reflect the effects of shock heating or radiative heating by a nearby 
star or protostar. As far as we are aware, all ISO observations of such sources (prior to those 
reported in this Letter) are consistent with an ortho- to para-H 2 ratio of three, the value expected 
at temperatures above ~ 200 K if the observed regions have been warm long enough to reach 
equilibrium between the ortho (odd J states with total nuclear spin 1) and para states (even J 
states with total nuclear spin zero){] 

In this Letter we report the results of ISO observations of molecular hydrogen in the source 
HH 54, a Herbig-Haro object located at the edge of the Chamaleon II dark cloud at an estimated 
distance ~ 200 pc from the Sun (Hughes & Hartigan 1992). The source of the outflow responsible 


1 Other earlier astronomical observations of H 2 also indicated ortho- to para-lR ratios consistent with the values 
expected in thermal equilibrium. Thus infrared absorption line observations of NGC 2024 IRS 2 (Lacy et al. 1994) 
implied an ortho- to para-H 2 ratio < 0.8, consistent with the value of 0.2 expected at the temperature (~ 45 K) 
of the absorbing gas, while ultraviolet absorption line observations toward hot stars (Savage et al. 1977) using the 
Copernicus satellite have revealed typical ortho- to para-LL ratios ~ 1 in diffuse clouds, again consistent with thermal 
equilibrium at the temperatures ~ 80 K that are typical of such clouds. We note that vibrational emissions from warm 
photodissociation regions and planetary nebulae often exhibit ortho- to para-Lh ratios in excited vibrational states 
that are smaller than the value of 3 expected in equilbrium at high temperature (e.g. Hasegawa et al. 1987, Ramsay 
et al. 1993, Chrysostomou et al. 1993, Hora & Latter 1996), with typical values in the range 1.7 - 2.1. However, 
Sternberg & Neufeld (1998) have argued that these ortho- to para-H 2 ratios measured for excited vibrational states 
do not require true ortho- to para-H 2 ratios smaller than 3 but are simply a consequence of optical depth effects in 
the fluorescent pumping of the observed vibrational emissions. 
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for HH 54 is believed to be IRAS 12496-7650 (Hughes et al. 1989), a deeply embedded young 
stellar object lying approximately 4' to the southwest. Previous observations of infrared H 2 
vibrational lines (Gredel 1994) and optical line emissions (Schwarz & Dopita 1980) have suggested 
the presence of hot, shock-excited gas in HH 54. 

Our observations are described in §2 below, and our results reported in §3 below. A discussion 
follows in §4, in which particular emphasis is placed upon the non-equilibrium ortho- to para-H 2 
ratio ~ 1.2 measured in HH 54. 


2. Observations and data reduction 

Using the SWS of ISO in its grating mode (SWS02), we observed the S(l), S(2), S(3), S(4) 
and S(5) lines of H 2 toward the source HH 54 on 1997 November 4th. The ISO beam was centered 
midway between two strong sources of H 2 vibrational emission, positions E and K (Sandell et al. 
1987), at coordinates a = 12h 55m 53.4s, 6 = —76° 56' 20.5" (J2000), with the long axis of the 
beam oriented at position angle 118°. The total observing time on target was 6400 s, including 
overheads for dark current measurements and calibration. The rest wavelength, spectral resolution, 
and beam size are given in Table 1 for each of the five observed H 2 lines. 


TABLE 1 

H 2 rotational lines observed toward HH54 E+K 


Line 

Wavelength 

Velocity 

Beam size 

Line Flux 6 

Line Intensity* 50 

Upper state 



Resolution 0 




column density 0 ^ 


Om) 

(km s -1 ) 

(arcsec) 

(W cm -2 ) 

(erg cm -2 s —1 sr —1 ) 

(cm -2 ) 

S(l) 

17.0348 

160 

14 x 27 

4.8 x 10~ 2 ° 

5.4 x 10“ 5 

1.3 x 10 19 

S(2) 

12.2786 

240 

14 x 27 

2.1 x 10 -19 

2.4 x 10“ 4 

7.5 x 10 18 

S(3) 

9.6649 

160 

14 x 20 

1.7 x 10' 19 

2.6 x 10~ 4 

2.2 x 10 18 

S(4) 

8.0251 

200 

14 x 20 

1.6 x 10" 19 

2.4 x 10~ 4 

5.0 x 10 17 

S(5) 

6.9095 

240 

14 x 20 

1.9 x 10” 19 

2.9 x 10~ 4 

2.2 x 10 17 


a FWHM of the instrumental profile for an extended source, from the ISO SWS Observer’s 
manual, Issue 3.0. For a point source, the spectral resolution is roughly a factor of two better. 
b uncorrected for extinction 
c beam averaged 

d computed with a correction for extinction (see text) 
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The initial data reduction was carried out with version 6.22 of the ISO pipeline software, and 
the ISAP software package]] was then used to remove bad data points and to co-add the individual 
spectral scans. 


3. Results 

Figure 1 shows the SWS H 2 spectra observed toward HH 54 E+K. In every case, the observed 
line width is no broader than the instrumental response function, i.e. the lines are unresolved at 
the resolution of the SWS grating mode. 

The measured line fluxes, beam-averaged line intensities and beam-averaged H 2 column 
densities in rotational states J = 3 through 7 are given in Table 1. The current SWS flux 
calibration is believed accurate to ~ 30%. In computing the H 2 column densities, we used the 
spontaneous radiative rates of Turner, Kirby-Docken, & Dalgarno (1977; confirmed recently by the 
more precise calculations of Wolniewicz, Simbotin & Dalgarno 1998) for quadrupole transitions 
of H 2 , which imply that the lines that we observed are all optically thin. We also applied dust 
extinction corrections based upon the interstellar extinction curves of Draine (1989), adopting 
Gredel’s (1994) estimate of 0.3 mag for E(J — H ) that was obtained from observations of the [Fell] 
lines at 1.257 and 1.644 /an toward HH 54 E+K. The assumed extinctions for the S(l), S(2), S(3), 
S(4) and S(5) lines were accordingly 0.10, 0.12, 0.30, 0.09 and 0.04 mag respectively]]]. 

Figure 2 shows a “rotational plot” for H 2 rotational states J = 3 through 7. Here the 
logarithm of Nj/(gjgs ) is plotted against Ej/k, where Nj is the beam-averaged column density, 
gj = 2J + 1 is the rotational degeneracy, gs is the spin degeneracy (1 for even J and 3 for odd 
J), and Ej is the energy of the state of rotational quantum number J. Open squares correspond 
to the values tabulated in Table 1, with column densities averaged over regions of size 14" x 27" 
for J = 3 and 4 and size 14" x 20" for J = 5, 6, and 7. The filled squares for J = 3 and 4 show 
values that have been corrected to the smaller beam size of 14" x 20", on the assumption that 
the emission region is small compared to either beam. That assumption, adopted hereafter in the 
discussion, is appropriate if the H 2 S(l) and S(2) emission follows the compact distribution of the 
higher-excitation vibrational emission mapped by Gredel (1994). 

Figure 2 shows that the points corresponding J = 3, 5, and 7 are colinear to within the 
observational uncertainties, consistent with the emission expected from gas in thermal equilibrium 


2 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and 
Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 

3 The assumed extinction for the S(3) line - and to a lesser extent for the S(1), S(2) and S(4) lines - is enhanced 
by silicate features in the assumed extinction curve (Draine 1989), the strength of which may vary from one source 
to another. However, we found that the gas temperature and ortho- to para-H 2 ratio derived below are only very 
weakly dependent upon the exact shape of the assumed extinction curve and are negligibly altered even if the silicate 
features are assumed to be entirely absent. 
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Fig. 1.— ISO Short Wavelength Spectrometer (SWS) spectra of pure rotational lines of H 2 observed 
toward HH 54 E+K, for a beam of size 14" x 27" [S(l) and S(2) lines] or 14" x 20" [S(3), S(4) 
and S(5)] centered at a = 12h 55m 53.4s, 5 = —76° 56' 20.5" (J2000). The S(2), S(3), S(4) and S(5) 
lines are respectively offset by 10, 20, 30, and 40 Jy. 
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Fig. 2.— Rotational plot for H 2 rotational states J = 3 through 7. The logarithm of Nj/{gjgs) is 
plotted against Ej/k, where Nj is the beam-averaged column density, gj = 2 J +1 is the rotational 
degeneracy, gs is the spin degeneracy (1 for even J and 3 for odd J) and Ej is the energy of the 
state of rotational quantum number J. Open squares correspond to the values tabulated in Table 
1, with column densities averaged over regions of size 14" x 27" for J = 3 and 4 and size 14" x 20 ,/ 
for J = 5, 6, and 7. Filled squares for J = 3 and 4 show values that have been corrected to the 
smaller beam size of 14" x 20", on the assumption that the emission region is small compared to 
either beam. Straight lines show the fit for a gas temperature of 650 K and an ortho- to para-H 2 
ratio of 1.2 (see text). 
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at a temperature T ~ 700 K. However, the column densities in J = 2 and 4 lie substantially above 
the best-fit line for the ortho (odd-J states), implying that the ortho- to para-H 2 ratio (OPR) is 
smaller than the high temperature equilibrium value of 3 expected at 700 K. Our best fit to all 
the data points implies an OPR of 1.2 ± 0.4, a gas temperature of 650 K, and a total H 2 column 
density of 9 x 10 19 cm -2 (averaged over a 14 /; x 20^ beam). 

The rotational temperature derived for HH 54 E+K is larger than the temperature ~ 330 K 
derived by Liseau et al. (1996) from ISO Long Wavelength Spectrometer (LWS) observations of 
CO transitions towards HH 54 B, a difference that presumably results from the different pointing 
and larger LWS beam size used by Liseau et al. On the other hand, the H 2 rotational temperature 
we derived towards HH 54 E+K is smaller than the temperature of 2100 K derived by Gredel 
(1994) from the observed ratio of the H 2 v = 1 — 0 and v = 2 — 1 vibrational bands. Again, such 
a discrepancy is unsurprising and probably indicates that a mixture of shock velocities and gas 
temperatures is present within the beam, with vibrational emissions selectively probing hotter gas 
associated with faster shocks. 


4. Discussion 

Our measurement of an OPR smaller than 3 in gas of temperature ~ 650 K is prima facie 
evidence that the gas has not been warm long enough to reach equilibrium between the ortho 
and para states; thus the measured OPR ~ 1.2 is the legacy of an earlier stage in the thermal 
history of the gas when the OPR had reached equilibrium at a value 1.2 (corresponding to a 
gas temperature ^ 90 K). 

We can place an upper limit on the time period, r, for which the emitting gas has been warm 
by estimating the timescale for conversion from para- to ortho-H 2 , r conv . In media of low fractional 
ionization, the conversion of para- to ortho-H 2 is dominated by reactive collisions with atomic 
hydrogen, for which the rate coefficient at 650 K is 1.0 x 10 _13 cm 3 s _1 (Tine et al. 1997): thus 
Tconv is given by 3000 [100 cm” 3 /n(H)] yr, where n(H) is the density of hydrogen atoms. (Even if 
other interconversion mechanisms are significant, this is in any case an upper limit on r conv .) 

Models for the chemistry of dense molecular gas (e.g. Neufeld, Lepp & Melnick 1995) 
predict that once the temperature of the gas exceeds ~ 300 K, atomic oxygen will be rapidly 
converted to water, a prediction that has been supported by recent observations of water toward 
warm shocked regions in the Orion Molecular Cloud (Harwit et al. 1998). The production of 
water takes place by means of a series of two hydrogen atom abstraction reactions with H 2 and 
therefore leads to the production of two hydrogen atoms for each water molecule produced; 
thus the atomic hydrogen density in recently-heated gas at temperature 650 K is at least twice 
the atomic oxygen density that was present in the gas prior to its being heated. Assuming 
an initial O abundance, n(0)/n(H2) of 3.5 x 10~ 4 (consistent with the gas-phase oxygen and 
carbon abundances of Cardelli et al. 1995 and the assumption that O accounts for most of the 



gas-phase oxygen not bound as CO), this places an upper limit of 5000 [10 4 5 cm 3 /n( H 2 )] yr on 
Tconv (This conservative upper limit still applies, of course, even if atomic hydrogen is produced 
by other mechanisms in addition to the reaction of O with H 2 to form water). Finally, the fact 
that the J = 7 state is apparently thermalized (c.f. Figure 2) implies a lower limit on the H 2 
density of ~ 10 5 cm ", the critical density at which the collisional de-excitation rate for J = 7 
(Tine et al. 1997) is equal to the spontaneous radiative decay rate (Turner et al. 1977). Thus 
we conclude that the emitting gas that we observed in HH 54 has been warm for a period 
t ~ T conv ^ 3000 [100cm _3 /n(H)] yr 5000 [10 5 cm 3 /n(H 2 )] yr ^ 5000 yr. 

Transient heating by shock waves - observed widely in Herbig-Haro objects - provides a 
natural explanation of the non-equilibrium OPR that we observed. Indeed, a non-equilibrium 
OPR within molecular gas heated by slow shocks has been predicted recently by Timmermann 
(1998), who presented a detailed theoretical treatment of ortho-para interconversion within slow 
‘C’-type shocks. Timmermann obtained predictions for the H 2 rotational line strengths for 
shocks of velocity 10 to SOkms^ 1 that propagate in gas of preshock H 2 density 5 x 10 3 , 5 x 10 1 
and 5 x 10 5 cm~ 3 and preshock ortho to para-H 2 ratio OPR* = 1 or 3. The results show that 
para-to-ortho conversion is incomplete for shock velocities smaller than 20 to 25 kms -1 - the 
timescale for equilibration being longer than the shock timescale — so that slow shocks with 
OPR* = 1 show a non-equilibrium OPR within the postshock region. 

Figure 3 shows a quantitative comparison between the observed S(2)/S(l) and S(3)/S(l) 
line flux ratios and Timmermann’s predictions for a preshock H 2 density of 5 x 10 5 cm -3 . In 
fact, none of the shock parameters actually considered by Timmermann yields a S(2)/S(l) ratio 
large enough (i.e. a postshock OPR small enough) to match the observations. However, we have 
interpolated and extrapolated Timmermann’s results to cases (dotted lines) where the initial 
OPR is other than 1 or 3, assuming that the line fluxes depend linearly^] upon the initial para-H 2 
fraction, (OPR* + 1) . Our extrapolation shows that a shock of velocity ~ 22kms -1 can provide 

a satisfactory fit to the S(2)/S(l) and S(3)/S(l) line ratios, given a sufficiently small initial ortho- 
to para-H 2 ratio, OPR,; 0.2, 

The results of Timmermann et al. (1998) also provide an explanation for why previous 
studies (e.g. Smith, Davis & Lioure 1997) of vibrational emissions from protostellar outflows and 
Herbig-Haro objects have always revealed an OPR close to 3: vibrational emissions inevitably 
trace faster, hotter shocks in which the atomic hydrogen fraction is enhanced by streaming 
ion-neutral collisions and the rate coefficients for reactive collisions between H and H 2 are larger. 
Thus Timmermann’s models predict that para-to-ortho conversion is rapid enough to yield an 
equilibrium OPR of 3 behind any shock that is fast enough to excite significant vibrational 
emission. 


4 A linear dependence would obtain exactly provided that the temperature structure of the shock is independent 

of the initial OPR. Our extrapolation method yields the expected result that the fluxes in ortho-IR transitions tend 
to zero in the limit of small initial OPR and small shock velocity. 
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Fig. 3.— Comparison between the observed S(2)/S(l) and S(3)/S(l) H 2 line flux ratios (filled 
square) and the predictions of Timmermann (1998) for a ‘C’-type shock propagating in gas of H 2 
density 5 x 10 5 cnr 3 , shown as a function of the shock velocity, v s , and the initial (i.e. preshock) 
ortho- to para-H 2 ratio, OPR*. The observed fluxes have been corrected for extinction (see text). 
We have extrapolated or interpolated Timmermann’s results to cases (dotted lines) where the initial 
OPR is other than 1 or 3, assuming that the line fluxes depend linearly upon the initial para-H 2 
fraction, (OPR* + l) -1 . 
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